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Reaction plane angle dependence of dihadron azimuthal correlations from a 
multiphase transport model calculation 
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Dihadron azimuthal angle correlations relative to the reaction plane have been investigated in An 
+ Au collisions at ^/snn = 200 GeV using a multi-phase transport model (AMPT). Such reaction 
plane azimuthal angle dependent correlations can shed light on path-length effect of energy loss of 
high transverse momentum particles propagating through the hot dense medium. The correlations 
vary with the trigger particle azimuthal angle with respect to the reaction plane direction, ^j, = 
<j)T — 'i'EP, which is consistent with the experimental observation by the STAR collaboration. The 
dihadron azimuthal angle correlation functions on the away side of the trigger particle present a 
distinct evolution from a single peak to a broad, possibly double peak, structure when the trigger 
particle direction goes from in-plane to out-of-plane of the reaction plane. The away-side angular 
correlation functions are asymmetric with respect to the back-to-back direction in some regions 
of <j)s, which could provide insight on testing vi method to reconstruct the reaction plane. In 
addition, both the root-mean-square width (Wrms) of the away-side correlation distribution and 
the splitting parameter D between the away-side double peaks increase slightly with 0s, and the 
average transverse momentum of the away-side associated hadrons shows a strong ^s dependence. 
Our results indicate that strong parton cascade and resultant energy loss could play an important 
role for the appearance of a double-peak structure in the dihadron azimuthal angular correlation 
function on the away side of the trigger particle. 

PACS numbers: 12.38.Mh, ll.lO.Wx, 25.75.Dw 



I. INTRODUCTION 



Lattice QCD calculations predicted that a novel state 
of matter distinct from ordinary hadronic matter, can 
be created in ultra-high energy heavy-ion coUisions. The 
BNL Relativistic Heavy Ion CoUider (RHIC) was con- 
structed to investigate the properties of the new state 
of matter . Studying the properties of this dense nu- 
clear matter is a great challenge both experimentaUy and 
theoreticaUy. Recently features of partonic collectivity 
and strongly interacting have been found for the mat- 
ter of extreme temperature and energy density created 
in nucleus- nucleus collisions at RHIC [j-Q- The energy 
loss of hard partons traversing the dense medium created 
in the collisions is a crucial question. Usually, there are 
two main mechanisms causing a jet to lose energy: elas- 
tic collisions with deconfined partons and induced gluon 
radiation 

[Ml, and the gluon radiation is considered 
as the main process responsible for energy loss. Both 
mechanisms suggest that the energy loss would depend 
strongly on the traversed path length of the propagat- 
ing jets [l3 - [l6| . Two particle azimuthal angular corre- 
lation has been demonstrated as a good probe to inves- 
tigate the interactions between jets and the hot dense 
medium [13 ■ It has been reported previously, at 
RHIC, that for a high pT trigger particle, the correlated 



yields at pT >2 GeV/c on the away side are strongly 
suppressed [l3|, while for lower px the yields are en- 
hanced. The correlated hadrons on the away-side appear 
to be partially equilibrated with the bulk medium and 
the azimuthal angular correlation distributions are much 
broader, possible with double peaks [20l - |22j . Recently 
many theoretical models have been proposed to describe 
the phenomenon: Mach-cone shock waves generated by 
large energy deposition in the medium (23j , conical emis- 
sion due to Cherenkov radiation the excitation of 
collective plasmon waves [Isj, lar ge angl e gluon radia- 
tion [H, [23] , strong parton cascade j28l - [3g | , jets deflected 
by radial flow [23, [Ml and path-length dependent en- 
ergy loss [13 ■ It has also been observed previously that 
the double peak structure in the dihadron angular cor- 
relation strongly depends on 0s, the trigger particle az- 
imuthal angle with respect to the direction of the event 
plane (33l . f34j . The angular correlation distributions on 
the away side are different between in-plane and out-of- 
plane orientations which indicate a strong path-length 
effect on parton- medium interactions (33 |. However, no 
Monte Carlo simulation has been performed to demon- 
strate such kind of reaction plane angle dependence so 
far. In this work, we will use a multi-phase transport 
model (AMPT) [33 to reveal the dependence of the reac- 
tion plane orientation on the dihadron correlation struc- 
ture which can be understood by the path-length effect 
of large elastic partonic energy loss. 
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AMPT model is a hybrid model which consists of four 
main components: the initial conditions, partonic inter- 
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actions, the conversion from partonic matter to hadronic 
matter and hadronic rescattering. AMPT model includes 
two versions: the default version and the string melt- 
ing version where strings are melted and hadronizations 
using parton recombination mechanism. Previous stud- 
ies have shown that the partonic effect could not be 
neglected and the string melting AMPT model which 
contains stronger parton cascade is much more appro- 
priate than the default AMPT model when the energy 
density is much higher than the critical density for the 
predicted phase transition [ssl - iSTj . The parton cascade 
model ZPC in the AMPT model only includes elastic 
2—7^2 partonic interactions [111, higher order inelastic pro- 
cess, which might become dominant at high densities dur- 
ing the early stage of relativistic heavy ion collisions [11] , 
have not been included. To fully model the parton en- 
ergy loss, in addition to the elastic energy loss the ra- 
diation ener gy l oss (many-body inelastic interactions in 
general) [Tol I 111] should also be implemented in this hy- 
brid transport model. However, it is not the intention 
of our current work to embark on such a comprehensive 
task of implementing radiation energy loss in AMPT. In 
order to match the parton energy loss in experimental 
data a large elastic scattering cross section (ct) of 10 mb 
has to be used in AMPT model. The large elastic cross 
section could contribute to the emergence of a double 
peak structure on the away side angular correlation dis- 
tribution. 

For mid-central collisions, the initial overlap region of 
the two nuclei looks like an almond. Figure 1 shows the 
initial parton coordinate distribution for 40% collision 
centrality from AMPT model simulations. The dimen- 
sions of dense matter distribution for in-plane and out- 
of-plane are quite different in mid-central collisions. We 
have calculated the eccentricity of the initial state, right 
after the collision in AMPT model, to be around 0.15. 
Particles along out-of-plane will on average propagate 
a longer path length than those along the in-plane di- 
rection. Triggered high momentum particles are biased 
preferentially from the surface area and outward from dif- 
ferent azimuthal angles. The corresponding away side jet 
and the associated particles will propagate through dif- 
ferent dense medium length. This can result in different 
energy loss due to dependence on interaction strength 
between the traversing parton and the dense medium, 
which may lead to different correlation structures on the 
back-to-back direction. 



II. ANALYSIS METHOD 

In order to investigate the path-length effect, we di- 
vide the whole initial angular phase space region into 16 
equal slices based on ~ (f>T — ep in Fig 2, where (pr 
and vPbp are the azimuthal angle of triggered particle 
and event plane, respectively. We calculate the dihadron 
azimuthal angle correlations between a high pT hadron 
(trigger hadron) and low pT one (associated hadron) sep- 



arately in each slice. The analysis method is similar to 
that used by the experiments [l^ I39l - l4l| . In the same 
event, pairs of an associated particle with a triggered par- 
ticle are accumulated to obtain A(j) = (f>assoc - i'trig raw 
distributions, where 4>assoc and (jitrig are the azimuthal 
angle of associated particle and triggered particle, respec- 
tively. In order to reproduce the background which is 
mainly from the anisotropic flow [13, 15], a mixing-event 
method is applied. We mixed the events which have very 
close centrality into a new event, and extracted Acf) dis- 
tribution to be used as a background distribution. A zero 
yield at minimum (ZYAM) assumption was adopted to 
subtract the background [41[. We can directly use the 
event plane angle ^f^p in our AMPT model simulation 
which sets ^ep = 0. 



III. RESULTS AND DISCUSSIONS 

Figure 3 shows the dihadron correlation functions for 
16 slices as determined by the trigger particle emission 
angle with respect to the reaction plane. The data are 
Au+Au coUisions at = 200 GcV from AMPT 

model simulation with a 20-60% collision centrality. We 
set pt range for trigger particles to be 2.5 < p^'^ < 6 
GeV/c and for associated particles to be 0.15 < p^^"" < 3 
GeV/c in our analysis. Both triggered and associated 
particles are further selected with a pseudo-rapidity cut 
of 1 77 1 < 1. We can observe distinct evolutions from in- 
plane to out-of-plane on both near and away side. The 
near side correlations maintain a single Gaussian struc- 
ture while the correlation yields decrease from in-plane to 
out-of-plane direction. On away-side, however, the cor- 
relation structure evolves from a single peak to double 
peaks. Furthermore, we note that the double peaks are 
not symmetric for data from some trigger particle emis- 
sion directions. For further discussion on the asymmetric 
structure, we focus on slice 2 and slice 15 as examples. 
We find that the asymmetry of their correlation struc- 
tures on the away side represents different orientation. 
For slice 2 as shown in Fig. 2, the height of peak cen- 
tered at around 2.8 rad is higher than that centered at 
around 3.8 rad, which is due to that on the away side 
direction particles travel through longer medium path 
length at around 3.8 rad (i.e. the left side medium of 
the back-side of slice 2, namely slice 10) than that at 
2.8 rad (i.e. the right side medium of the slice 10). For 
slice 15, it shows the opposite asymmetry. To guide the 
eyes, the two Gaussian fits are plotted for the Slice 2 
and 15. Consequently, hadrons travel through different 
medium length and induce different asymmetry of the 
away-side dihadron structure. In other words, the cor- 
relation structure depends on the path-length that the 
associated particles travel. 

Our simulation result seems to indicate that violent 
parton cascade processes can make significant contribu- 
tions to the splitting structure in the away-side dihadron 
correlations. Recent preliminary results by Ma et al. 
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showed that the double peak structure of final-state di- 
hadron correlation can also be generated when a back jet 
traverses the dense core due to the shadowing effect of 
the core [i^. In this picture, the dense core shadowing 
effect is induced by the transverse expansion or radial 
flow of colliding system. When partons traverse through 
a dense core region, they will suffer stronge interactions 
with the medium, which will lead to extensive diffusions 
and particle redistributions. In the event-averaged case 
when jet traverses the dense medium, the double peak 
structure could be due to the sum of deflected and/or 
cone-like jets. This dense core shadowing effect is also 
related to path length dependence of jet-medium inter- 
actions. The longer path length direction coincides with 
the direction into the dense core, jet particles are blocked 
strongly due to the shadowing effect. The asymmetry in 
the path length distribution for particles interacting with 
the dense core medium will result in the asymmetry in 
the dihadron correlation. 
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FIG. 2: (Color online) Schematic illustration of 16 slices of 
trigger particle azimuthal angle relative to the event plane, 

(j>s = (t>T - EP- 
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FIG. 1: (Color online) The initial parton distribution in trans- 
verse plane from AMPT model simulation of Au+Au colli- 
sions at y/SNN = 200 GeV for 40% centrality. 

To further investigate the path-length effect, we com- 
pared correlation functions from mid-central collisions 
with central collision results. We show dihadron az- 
imuthal angle correlation functions for 0-10% centrality 
and 20-60% centrality from AMPT model simulations in 
Figure 4. We divide the flrst quadrant of initial phase 
space region into 6 equal slices based on 0^ and in order 
to increase the statistics, we fold the other three quad- 
rants into the first quadrant. Significant differences in the 
correlation functions between these two centrality bins 
become apparent. For mid-central collisions, we observe 
an evolution from a single peak to double peak struc- 
ture. But for central collisions, the correlation functions 



show little variation with 4>s and double peak structure 
is present for all slices. This can be explained that for 
central collisions, the overlap region of the two collid- 
ing nuclei is almost spherical and the dimensions of the 
dense medium are such that both in-plane and out-of- 
planc path length are long enough to generate a splitting 
structure. We conclude that in order to generate a double 
peak structure, a dense medium of sufficiently large size 
is needed [1^. Our simulation results indicate that the 
onset of the splitting structure corresponds to (ps ^ 45° 
for 20-60% centrality. The associated particle yields per 
triggered particle for central collisions is larger than that 
for mid-central collisions on both near and away side. 

In what follows, we will extract a few parameters from 
the dihadron angular correlation functions to describe the 
away-side correlation shape and amplitude. Systematic 
comparisons of these parameters between experimental 
results and theoretical calculations will provide insight 
on the response of the medium to the jet energy loss, and 
may provide new constraints for understanding compet- 
ing mechanisms for energy transport. 

Figure 5 presents the AMPT model calculations of the 
trigger direction dependence of the away side broaden- 
ing from Au-|-Au collisions at ^/snn = 200 GeV. For 
central collisions, the away-side width Wrms grows very 
slightly with (pg. For 20-60% central collisions, however, 
the away-side Wrms grows distinctly with trigger direc- 
tion (j)s which qualitatively agrees well with the experi- 
mental measurement [s^l ■ For (ps close to rad, the Wrms 
for 0-10% Au-|-Au collisions is already large enough which 
represents a remarkable broadening of the away-side cor- 
relation function, while that for the 20-60% Au4-Au col- 
lisions is relatively small corresponding to a narrow cor- 
relation. This is consistent with expected dependence on 
the traversing path length in the reaction plane direction 
between these two centrality bins. For (ps close to 7r/2 
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FIG. 3: (Color online) Background subtracted dihadron cor- 
relations with trigger particles in 16 slices of azimuthal angle 
relative to the event plane, 4>s = (f>T - ep- The data are from 
the melting AMPT model simulations with hadron rescatter- 
ing on Au+Au collisions at ^/JWn = 200 GeV for 20-60% 
centrality. (pT window cut in model: 2.5 < p^'® < 6 GeV/c 
and 0.15 < p^^"" < 3 GeV/c). Both the trigger and associ- 
ated particles are restricted within I77I < 1. The solid line in 
Slice 2 and 15 represents the two-Gaussians fit for the away 
side dihadron correlation. 

rad, the Wrms of these two centrality bins are very much 
close to each other. This can be attributed to the fact 
that the path length perpendicular to the reaction plane 
is not much different between these two centrality bins. 
The path length variation with the trigger particle ori- 
entation in the central collisions is relatively small com- 
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FIG. 4: (Color online) Background subtracted dihadron cor- 
relations for trigger particles of 2.5 < p^*® < 6 GeV/c and 
associated particles of 0.15 < Vt"° < 3 GeV/c from AMPT 
model simulation of Au+Au collisions at ^snn = 200 GeV 
with trigger particles in 6 equal slices of azimuthal angle from 
the event plane in the first quadrant for 0-10% collision cen- 
trality (open circles) and 20-60% collision centrality (solid 
circles). The rj window is restricted within |?7| < 1. Lines 
represent the two-Gaussians or one-Gaussian fit for away-side 
correlation. 



pared to that in the mid-central collisions. Accordingly, 
the away side Wrms for central collisions varies little with 
(j)s and is larger than that for mid-central collisions at all 
trigger particle orientations. 

Figure 6 shows the trigger particle azimuthal angle de- 
pendence of the away-side splitting parameter D (half 
distance between two splitting peaks on away-side in di- 
hadron azimuthal angle correlation). For mid-central 
collisions, the splitting parameter D grows significantly 
with the trigger particle azimuthal angle while for cen- 
tral collisions it grows only slightly. It can also be at- 
tributed to different path-length variations between these 
two centrality bins. We also show the trigger particle az- 
imuthal angle dependence of correlation yields for both 
near and away side in Figure 7. On near-side, the yields 
decrease with trigger azimuthal angle while on away-side, 
the yields increase. 

Figure 8 presents the A0 dependence of mean pT for 
in-plane and out-of-plane associate particles. We have 
combined data of the slices 1, 2, 7, 8, 9, 10, 15 and 16 into 
in-plane result and the rest into out-of-plane result. On 
the near-side, the mean pT values for both in-plane and 
out-of-plane agree well with each other and they stay al- 
most constant. On the away-side, they show different an- 
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FIG. 5: (Color online) The away-side Wrms of the diliadron 
correlation function versus the trigger particle azimuthal an- 
gle relative to the event plane, 03, from the AMPT model 
calculation with pt window cut (2.5 < p^** < 6 GeV/c and 
1 < Pt""" < 1.5 GeV/c ) in 0-10% (open circles) and 20-60% 
(solid circles) Au-I-Au collisions at ^sjvjv =200 GeV. The rj 
window is restricted within [r;! < 1. The preliminary experi- 
mental trend is for 20-60% centrality [3^ . 
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FIG. 6; (Color online) The splitting parameter 7? as a func- 
tion of trigger particle azimuthal angle from the event plane 
for trigger hadrons of 2.5 < p^** < 6 GeV/c and associated 
hadrons of 1 < p^r'""" < 1.5 GeV/c from AMPT model sim- 
ulations for Au+Au 0-10% (open circles) and 20-60% (solid 
circles) collisions at ^ysNN ~ 200 GeV. 

gular dependence. The meanp^ for in-plane stays almost 
flat while for out-of-plane, the mean px displays eoncave 
structure with the minimal value at around tt rad. Pre- 
viously published results indicate that for central colli- 
sions the concave structure implies that harder associated 
hadrons prefer larger angles with respect to the back-to- 
back direction and it may be attributed to the average 
interaction length of the away-side partons interacting 
with the medium is maximal at A(f> = tt rad (4^ . Our 
results indicate that for mid-central collisions, the path- 



length on away-side for in-plane is relatively small, and 
it varies very little with A0. For out-of-plane, however, 
partons on the away-side direction will travel through a 
long enough path length in the medium leading to a sig- 
nificant energy loss; partons emitted away from the back- 
to-back direction will traverse through relatively smaller 
path length with a correspondingly smaller energy loss. 
The path-length effect plays an important role on the 
structure of mean pT for the away side particles. 
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FIG. 7: (Color online) The yields of the dihadron corre- 
lation function on near (solid circles) and away side (open 
circles) versus the trigger particle azimuthal angle from the 
event plane, ^s, with pr window cut (2.5 < '® < 6 GeV/c 
and 1 < p^'"'"" < 1.5 GeV/c) in 20-60% Au-hAu collisions at 
^/sNN =200 GeV. The 77 window is restricted within |»7j < 1. 
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FIG. 8: (Color online) The mean ""^ versus A(j> {(fjassoc — 
4>trig) for in- (open circles) and out-of-plane (solid circles) in 
dihadron correlations with pT window cut (2.5 < p^ < 6 
GeV/c and 0.15 < pr""" < 3 GeV/c) in Au+Au collisions 
(20-60%) at ^/slTiv =200 GeV. The rj window is restricted 
within 1 77! < 1. 
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IV. CONCLUSIONS 

In conclusion, dihadron azimuthal angle correlation 
functions with respect to reaction plane have been stud- 
ied in the framework of a hybrid dynamical transport 
model, namely the AMPT model. In this model, the 
energy loss is essentially from elastic parton-parton scat- 
tering process, and the gluon radiation energy loss mech- 
anism has not been included yet. In order to match the 
experimental observations, a large elastic cross section 
has to be used. With the AMPT model simulations we 
observed the emerging of the splitting structure in the 
away-side dihadron correlations from Au+Au collisions 
at 200GeV/c. Wc found that the structure in the az- 
imuthal angle correlation functions has a strong depen- 
dence on the trigger particle orientation relative to the 
reaction plane, which is consistent with the observation 
from the STAR experiment. The correlation amplitude 
on the near-side slightly decrease from the in-plane to 
the out-of-plane direction, while the shapes of the away- 
side angular correlation functions evolve from a single 
peak to double peaks. For the triggered particles emit- 
ting from middle azimuthal angles relative to the reaction 
plane, the corresponding away-side correlation function 
is asymmetric which indicates that the shape of the corre- 
lation shall be sensitive to the traversing path length. By 
comparing results in the top 0-10% and 20-60% collision 
centrality bins, we found that the away-side path length 
in the reaction plane direction in 20-60% Au-|-Au collision 
was not long enough to generate a double peak structure 
while that in the top 0-10% collision it was long enough 
to do that. While in the out-of-plane direction, the cor- 
relation structure seems to be not very different since the 
path length is similar between these two centrality bins 
along this direction. Wc extracted a few parameters from 
correlation functions and discussed the possible underly- 
ing mechanisms for the splitting structure. Our results 
indicate that the away-side Wrms and splitting parameter 
D increase with (j)s for 20-60% centrality while for 0-10% 
centrality they change very slightly. This observation un- 
derscores the importance of the path-length effect when 



the away-side associated particles traverse the medium. 
The yields on near-side decrease slightly with 0s while 
on the away-side they increase. We also investigate the 
A(f> dependence of associated particles mean px- The 
dip structure of mean px as a function of Acj) along the 
out-of-plane direction is more prominent than that along 
the in-plane direction. These results could shed light on 
path-length effect of clastic collision energy loss for high 
transverse momentum particles propagating through the 
hot and dense medium. Strong parton cascade processes 
and resultant energy loss play important roles in the in- 
teraction between the jet-like particles and medium in 
the AMPT model. Furthermore, the preferential way to 
reconstruct the reaction plane in experiment now is using 
the m-th harmonic plane, of which the second harmonic 
plane, i.e., ?j2-mcthod. is often used (s^ . If one uses 
the first harmonic plane {vi method), it will be possi- 
ble to distinguish between structures of different trigger 
particle emission regions, thus allowing to separate dif- 
ferent path length effect in experiment. We note that the 
present findings are mostly relevant for the elastic energy 
loss and we do not know contributions from gluon radia- 
tion loss and their consequence on dihadron correlations 
yet in a dynamical model simulation. Nevertheless, our 
studies indicate that path length effect in jet medium in- 
teractions is very important and should be thoroughly 
investigated both experimentally and theoretically. 
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